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The unsolved problem of determining which densities are ground state densities of an interacting 
electron system in some external potential is important to the foundations of density functional 
theory. A coarse-grained version of this ensemble V-representability problem is shown to be thor- 
oughly tractable. Averaging the density of an interacting electron system over the cells of a regular 
partition of space produces a coarse-grained density. It is proved that every strictly positive coarse- 
grained density is coarse-grained ensemble V-representable: there is a unique potential, constant 
over each cell of the partition, which has a ground state with the prescribed coarse-grained density. 
For a system confined to a box, the (coarse-grained) Lieb functional is also shown to be Gateaux 
different iable. All results extend to open systems. 



1. INTRODUCTION 

Four decades ago, Hohenberg and Kohn launched 
modern density functional theory with their famous 
theorem [l| stating that a density function for a sys- 
tem of interacting particles can be a ground state den- 
sity of at most one external potential. A density is 
said to be ensemble V-representable (EV-representable) 
if there is a mixed ground state (i.e., a density ma- 
trix) for some potential which has that density. The V- 
representability problem, of determining which densities 
are EV-representable not only remains unresolved to this 
day, but has seen almost no progress. 

Optimistically, one might have hoped that any density 
which can arise from a state with finite energy is EV- 
representable. But there are known counterexamples [H, 
0] based on nonanalyticity. It is also easy to con- 
struct densities, by modifying the short wavelength parts, 
which can come only from states with infinite kinetic en- 
ergy, and thus certainly cannot be ground state densi- 
ties. These observations point to the potential useful- 
ness of suppressing short length-scale degrees of free- 
dom, or otherwise rendering them innocuous. Formu- 
lating the theory on a lattice achieves that by eliminat- 
ing short wavelengths. The lattice program was initi- 
ated by Kohn0| and Chayes, Chayes and Ruskai[H later 
proved that all densities on a lattice are ensemble V- 
representable@, 0, 0| ■ 

Unfortunately, working on a lattice fundamentally al- 
ters the problem, and it is difficult to relate the results to 
standard continuum quantum mechanics. In this paper, 
I propose a way to make the short length-scale degrees of 
freedom harmless without altering the underlying quan- 
tum mechanics. Instead, by coarse- graining, we impose 
a bound on the precision with which we permit ourselves 
to specify a desired density. Specifically, we divide our 
space, either all of M'' or a parallelepiped, into a regular 
array of parellelepiped cells of volume fl. The cells are 
labelled Si, i?2, • ■ • (a countable infinity), or Bi, . . . , Bm, 
according to the case. Denoting the operator correspond- 
ing to the number of particles in cell Bk by Nk, we get 



a coarse-grained density operator p, with components 

Pk = Nk/n. 

Now we can ask about states of our system of parti- 
cles which have a specified coarse-grained density p. In 
particular, we want to know whether there is a potential 
which has a (mixed) ground state with coarse-grained 
density p. If so, p is said to be coarse-grained EV- 
representable (cg-E V-representable). As will be shown, 
all coarse-grained densities, except those which put zero 
particles in one or more cells, are coarse-grained EV- 
representable. Moreover, among the potentials which are 
constant on each cell, there is one and only one having 
a ground-state coarse-grained density p. The proof of 
coarse-grained EV-representability is the central result 
of this paper. It should be clear, but bears emphasizing, 
that among all the states yielding the coarse-grained den- 
sity p, the one which is guaranteed to be V-representable 
by this result minimizes the kinetic-plus-interaction en- 
ergy. In the case of a system of particles confined to 
a box, we show that Fj^, viewed as a function on the 
affine hull of the set where it is finite, is Gateaux dif- 
ferentiable. In what sense can be differentiable at 
EV-representable densities is a question of continuing 
interest In the unconfined case, Fl is not Gateaux 
differentiable anywhere, at least in the ordinary sense. 

Coarse-graining allows short length-scale degrees of 
freedom to do as they like, without artifically remov- 
ing them. Arguably, it is also more harmonious with 
practical computational limitations than the traditional 
approach. 

In the standard (fine-grained) theoryQ, it is known 
that the densities which are EV-representable are dense 
in the effective domain of Fl (see Thm. 3.13 of Lieb@ 
or Sec. 1.6 of Ekeland and TemamQ). The useful- 
ness of this observation is limited by the absence of any 
criterion to pick out the good densities, and the fact 
(relevant for Kohn-Sham theory[l3|) that they may de- 
pend on the particle interaction. Within its scope, the 
coarse-grained theory does not suffer from these prob- 
lems. One might suspect that the coarse-grained EV- 
representability statements follow trivially from the ex- 
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istence of that dense set of densities, but it is not so be- 
cause the set of fine-grained densities belonging to a given 
coarse-grained density does not contain any neighbor- 
hood. 

The body of this paper unfolds as follows. Section [2] 
describes, in a rough way, the main ideas involved and 
the problems to be grappled with. Details of the confined 
case (in a box) are given in Section [3l Technically, this 
case is much easier than the unconfined case. Thus, al- 
though it stands on its own, it also serves as a warm-up 
exercise for the unconfined case, which is treated in Sec- 
tion [H In Section [U the extension to open systems with 
indefinite total particle number is discussed, in somewhat 
less detail. All the results carry over, and the requisite 
modifications of the proofs are sketched. Finally, in sec- 
tion [6l the results are put into perspective and directions 
for future work indicated. 

In the body of this paper, the word density without 
modifier will carry the usual meaning, namely, "fine- 
grained density." Frequently, coarse-grained will be ab- 
breviated as eg. Coarse-grained density and number are 
usually expectation values; when used in their operator 
aspect, a caret will be added {Ni etc.) to help avoid con- 
fusion. This practice is not extended to other operators, 
for which ambiguity is not a danger. 



2. MAIN IDEAS 

This section is an impressionistic sketch of the main 
ideas involved in the technical work of the following three 
sections, concentrating on the unconfined case with fixed 
particle number. As described in the introduction, we 
cover M.^ with identical cells. All (fine-grained) densi- 
ties with the same average density in each cell {pi in Bi) 
are gathered together into an equivalence class called a 
coarse-grained density. A cg-density p = {pi, p2, . . .) sat- 
isfies Pi > and J2i Pi — ^ ■ Dual to the cg-densities, 
we introduce a space V of external potentials which are 
constant on each cell; other restrictions will be imposed 
later. For such a potential, the integral J vp dx depends 
only on the coarse-grained class of the density p. 

The Lieb functional is a functional of coarse-grained 
density which represents the miminal intrinsic energy 
(kinetic plus Coulomb interaction) consistent with the 
specified cg-density: 

^^l(p) =inf{Trrffo : T ^ p} , 

where Hq is the kinetic-plus-interaction Hamiltonian and 
the customary notation F ^ p indicates that mixed state 
F has cg-density p. Fj^ is not defined as a minimum be- 
cause it is not obvious that there is any state which re- 
alizes the infimum; finite-dimensional intuition is a bad 
guide here. But only a state with minimal intrinsic en- 
ergy in its class can possibly be a ground state for a po- 
tential from V, so it is crucial for cg-EV-representability 
that we show F]^ really is a minimum. 



To illustrate the kind of problem which can arise 
in ensuring the existence of a minimum in infinite- 
dimensional situations, consider trying to minimize the 
function f{-ip) = X]'^i(l/")l(V'l'^n)P over the unit vec- 
tors in a Hilbert space, where {</>„} is an orthonormal 
basis. No matter how small e > 0, = {/ < e} is 
infinite-dimensional (and yet ClSe is empty). Absent de- 
tailed knowledge about the intrinsic energy, the way to 
show that the Levy search^ implicit in the definition 
of i^L has a solution is to show that the search is "al- 
most finite-dimensional." For Fl, the gist of the argu- 
ment may be phrased semiclassically: Rough knowledge 
about localization of the density approximately confines 
the search to a bounded region of configuration space, 
while considering only states with relatively low intrin- 
sic energy gives localization in momentum. The two to- 
gether approximately confine the search to a bounded 
volume of phase space. Since a volume (27r?i)^^ of phase 
space corresponds to one dimension in Hilbert space, the 
search is almost-finite-dimensional. 

The ground state energy E{v) of v £ V is the lower 
bound on energies of states in the presence of v, which 
can be expressed as 

E{v) = inf |fl(p) + J vpdx^ . 

(Note that, despite the name, E{v) may fail to be at- 
tained, as for the potential which is everywhere zero.) 
This equation shows that E is the Legendre transform of 
Fl. It also shows that, for given p, 

F{p) + J vpdx> E{v) (1) 

holds for all v. Now, if the cg-density p is cg-EV- 
representable, then this inequality must be an equality 
for the realizing potential, and conversely, if it is equality 
for some v, then p is clearly cg-EV-representable. 

So, the task is to show that equality is achieved in in- 
equality ID). To that end, more structure is useful. We 
embed the cg-densities in the vector space X = £^ oi real 
sequences p such that = ^ \pi\ < oo, which is a Ba- 
nach space equipped with the norm ||p|| . Then we declare 
that we consider only potentials in X*, the topological 
dual of X, which consists of v which are bounded. The 
double Legendre transform of -FL is its lower semicontin- 
uous convex envelope. Therefore, by showing that -FL is 
lower semicontinuous and convex, we establish that Fl 
is (up to signs) the Legendre transform of E: 

F{p) = sup |ii'(i') — J vpdx}j . 

Convexity is immediate from the definition, but lower 
semicontinuity requires some work. 

In the case of a system confined to a box, X and X* 
are finite-dimensional. In that case, existence of a po- 
tential which gives equality in Eq. ([T]) is an automatic 
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consequence of the mutual Legendre transform relation 
oiE and Fl. 

The unconfined case is harder. We attack it by noting 
that the second Legendre transform relation says that 
equality in Eq. ([T]) can be approached as closely as de- 
sired. That is, a sequence of t> G X* can be found such 
that 

E{Vn) - j Vnpdx F{p). 

The trick is then to show that this sequence Vn converges 
in some sense and that the limit v^o satisfies equality 
in ll]). In carrying that out, it is necessary, in gen- 
eral, to locate the limit outside of X* entirely. This 
may seem surprising, at first, but was entirely to be ex- 
pected. After all, the harmonic oscillator ground state is 
v-representable, but the harmonic-oscillator potential is 
not bounded, hence not in X*. 

This, in brief, is how cg-EV-representability is es- 
tablished. We also prove a coarse-grained version of 
the Hohenberg-Kohn Theorem, with only a little more 
work than for the usual version. With that, we can 
see that, not only is every everywhere-positive p cg- 
EV-representable, but that the representing potential is 
unique. 

3. CONFINED SYSTEMS 

In this section, we treat confined systems. Standard 
tools of convex analysis will be used, without much com- 
ment apart from reminders of definitions. The reader for 
whom the material is unfamiliar can find it in a num- 
ber of textbooks, for ex amp le, chapter I of Ekeland and 
TemamQ, or van Tiel[TT|. Rockafellarfl^ treats the 
finite-dimensional theory exhaustively. In citing specific 
results in those works, the abbreviations ET, vT and R 
are used. 

We are concerned with a system of N fermions ( "elec- 
trons," with or without spin) in a box, which we take to 
be a parallelepiped for simplicity. The Hilbert space of 
antisymmetrized states for these particles is denoted TL. 
A mixed state, not necessarily normalized, has the form 
r — X^J^i "^«lV'n)('0n|, where ^„ are normalized states 
in 7i, < a„ for each n, and Yl°^=i < oo. The set 
of normalized mixed states, which are those satisfying 
TrF = X^J^i*^™ — 1; ^ill t)e denoted simply by S. The 
mixed states belong to the space Ci{TC) of trace class 
operators on TC. More specifically S is the intersection 
of the convex cone C^iji.) of positive operators with the 
hyperplane TrF = 1. Ci{H) is a Banach space under 
the trace norm ||A||fr — Tr|A|, where \A\ is the unique 
positive square root of A^A. 

The box is subdivided into a regular array of paral- 
lelepiped cells Bi, - ■ ■ , Bm of volume ft, with associated 
density operators Ni/il, ■ ■ ■ ,N]\j/il. Generally, we look 
upon sequences of this sort as points in K*^, and denote 
them by bold-face symbols. S decomposes into strata 



labeled by p £ R^^: 

S{p) = {r eS: TrTiV = p}. (2) 

The space into which cg-densities have been embedded is 
generically denoted by X. In this section, X — R*^. The 
term eg- density is reserved for elements of X which have 
only nonnegative components, denoted by X^. The set 
of normalized cg-densities is 

X+ ^i^peR^' : p>0,J2p,n^Ny (3) 

If p is not a normalized cg-density, then S{p) is empty. 
X^ is an (M — l)-dimensional simplex contained in the 
M ~ 1 dimensional hyperplane Xn of p G R^^ such that 
Pi — The interior of X^ relative to Xn corre- 
sponds to normalized cg-densities which are greater than 
zero in every cell, and will be denoted by X^~^ . (Stan- 
dard nomenclature [T2| calls Xn the affine hull of X^, 
and X^'^ the relative interior oi X^ .) 

The Hamiltonian Hq includes kinetic energy and the 
interaction between the electrons, with periodic, Dirich- 
let or Neumann boundary conditions. We assume the 
interaction is the repulsive Coulomb interaction, possi- 
bly modified by a bounded interaction, such as would 
be required to adapt it to periodic boundary conditions. 
The bounded addition may even be a multiparticle in- 
teraction, but we assume it has a constant adjusted so 
that the total interaction is nonnegative. i/o is then un- 
bounded, but bounded below by zero. 

There is a lower limit to the kinetic plus interaction 
energy of states having a given cg-density, the Lieb func- 
tional 

:= inf {TrTi/o : T G S{p)] . (4) 

If S(p) is empty, Fi^{p) — oo. It may appear at first sight 
that the trace in the definition is potentially not well- 
defined, since Hq is unbounded. However, since Hq > 
and r is a positive operator, there is no real difficulty. We 
can, for example, replace Hq by HqPhq (A), where Phq (A) 
is the spectral projection onto Hq < A, and take the limit 
A ^ oo. The set of states which realize the infimum in 
Eq. (g]) is denoted S{p): 

S{p) := {r G S{p) : TrEiJo = F^{p)} . (5) 

We will show later that this set is nonempty. Fl is 
bounded both above and below: < Fl < Fl""^. The 
lower bound is simple, and the upper bound follows by 
using the fact (A^-representability) that a smooth density 
with any desired cg-density can be realized by an appro- 
priate mixed state, and then using the bounds on Fl in 
terms of Vp as derived in LiebQ (see Thms. 3.8 and 3.9 
there) . The largest of these (finite) bounds is for the case 
of all particles in a single cell, so Fl is bounded above. 

Having dealt with densities, we now turn to potentials. 
The potentials we want to consider take the constant 
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value Vi on cell Bi, information which is encoded in the 
vector V = {vi, . . . , vm) G K*^. We will write V for this 
space of potentials. Adding the potential v to Hq gives 

M 

H{v) = Ho + Y^ flv.p,. (6) 

When convenient, we use the alternate, and suggestive, 
notations / v p{x) dx ot v ■ p for the potential term. 
The ground state energy for potential G V is 

E{v) ■.= ini{TT:TH{v) : T e 5} . 

This really is the energy of some state; the infimum is 
sure to be realized in the present setting, as follows from 
the fact that our system is in a box (cf. Thm. 13. 2p . We 
can rewrite E(y) in terms of the Levy constrained search 
as 

E{v) = mi{Fi^{p) + V ■ p} (7) 
p 

That is, —E{—v) is the Legendre-Fenchel transform (also 
known as convex conjugate^ conjugate, or polar) of Fi^(p). 
Thus, E{v) is concave and upper semicontinuous where 
it is finite. But, the infimum in Eq. ([T]) is certain to lie 
between min{wi} and F™°-^ + ma.x{vi}, so that E{v) is 
finite everywhere on R*^, and therefore continuous. The 
Legendre-Fenchel transform of —E{—v) is the closed con- 
vex hull of i^L- is easily seen to be convex from its 
definition, so continuity on X^'^ is automatic because of 
finite dimensionality (ET: Cor. 1.2.3, vT: Thm. 5.23, 
R: Thm. 10.1), but lower semicontinuity at the relative 
boundary of is not. In fact, _Fl is lower semicontin- 
uous there as well. Rather than provide the proof (see 
Thm. 14. 3p , we will here simply assume that .Fl is defined 
on the relative boundary so that it is lower semicontinu- 
ous. This is harmless because the cg-densities there are 
not of interest. We summarize the foregoing discussion. 

Proposition 3.1 (convex duality). The functions F\^{p) 
and E{v) on R*^ are in duality, i.e., they are (up to 
signs ) Legendre transforms of each other. Hence, Fi^ and 
—E are convex and lower semicontinuous. 

One aspect of this relationship between Fj^ and E 
which will be important later is that G V is regarded 
as a linear functional on R*^ via p ^ Jvpdx. Since 
R^^ is finite-dimensional, v is automatically a continu- 
ous functional (though v{x) is not a continuous function 
of x unless it is constant). 

Now, we are ready for a coarse-grained version of the 
Hohenberg-Kohn theorem. If potentials v and v' differ 
by an overall constant, write ?; ~ d', and denote the 
equivalence class of v in V /R by [v] . 

Theorem 3.1 (eg Hohenberg-Kohn). Given S{p), 
there is at most one equivalence class in V/R which has 
a ground state in S (p) . 



Proof. First, v ^ v' implies that Fl(v) and H{v') have 
disjoint ground state manifolds in Ti. This is a result 
about partial differential equations and is proven just 
as for the conventional Hohenberg-Kohn theorem, being 
just a specialization to potentials in V. 

If G V has a ground state in S{p), then that state is 
in S_{p), and all elements of S_{p) are ground states for 
V, since all such states have the same value of TriJoF. 
Thus, if both v and v' have ground states in S{p), then 
they share a ground state, which contradicts the previous 
paragraph. □ 

This proof carries over to the unconfined setting of 
Section [5] and the Fock space setting of Section O with 
no essential changes. 

As just observed, the only states in S{p) which can 
possibly be ground states are those in S{p). The next 
task is therefore to show that S_{p) is nonempty. 

Theorem 3.2 (Existence of Levy search solutions). 

For every p G 

'^N' ^(P) / 0- That is, there is a normal- 
ized mixed state F such that F{p) = Ti HqT. 

We prove this by means of a couple of lemmas. It 
will be useful to write the intrinsic energy in a functional 
notation: 

f(F) :=TrFiJo, (8) 

regarded as a function on S, so that £^^{0,M] consists 
of normalized mixed states with intrinsic energy not ex- 
ceeding M. 

Lemma 3.1. £ is lower semicontinuous with respect to 
trace norm on S. 

Proof Suppose F„ G £-^{0,M] for n = 1,2,... and 
F„— >r in trace norm. We need to show that f (F) < M. 
Suppose instead that £(F) > M. 

If Pffo(e) denotes the spectral projection onto Hq < e, 
Tr THqPh„ (e) is monotonic in e so that there is some fi- 
nite e for which Tr FiJo-Pffo ("s) > But, trace-norm 
convergence implies weak convergence and HqPhq (e) is 
a bounded operator. Therefore, Tr TnHoPHo (e) exceeds 
M for large enough n since it converges to Tr THqPho (e). 
This, however, is impossible because M > 5(F„) > 
TrF„i/oPffo(e). □ 

This result is very general, depending as it does only on 
the positivity of mixed states and the semiboundedness 
of the Hamiltonian. The theorem and proof apply to the 
unconfined case, and even in Fock space. 

For the next lemma, we recall the concept of total 
houndedness. A set in a metric space is called totally 
bounded if, given e > 0, it can be covered by a finite 
number of balls of radius e. For a subset of a complete 
metric space, total boundedness is equivalent to relative 
compactness, i.e., having a compact closure. A Banach 
space is a complete metric space. 

Lemma 3.2. £~^(0, M] for M < oo is relatively compact 
with respect to trace norm. 
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Proof. £~^{0,M] is certainly bounded since it is con- 
tained in the unit ball of C^{H). The strategy is to 
show that, given e, there is a finite dimensional space W, 
such that all points of £~^{0,M] are within distance e 
of W. Since a bounded subset of W is totally bounded 
and e is arbitrary, this will show that £^^{0, M] is totally 
bounded, hence relatively compact. 

The interaction energy is nonnegative, so £~^{0,M] 
can only get bigger if we drop the interaction. Thus, it 
suffices to prove the lemma for the case that Hq ~ T 
contains only kinetic energy. 

Due to our assumptions, the single-particle eigenstates 
ipm of T can be written down explicitly, and the spectrum 
is certainly discrete. Up to energy M/e, there are only a 
finite number, tpi, . . . , tpx, and the set of Slater determi- 
nants made from these likewise span a finite-dimensional 
space, W. 

Now, any state T can be written as a linear combina- 
tion of Slater determinants of the (pm- But, if £{T) < M, 
the terms containing (pm for m > K must, all together, 
have norm less than e. That is, F is the sum of something 
in W and something with norm less than e. □ 

The box played a crucial role in this proof; its counter- 
part in the unconfined case will be much more difficult. 

Proof, (of Theorem I3.2p Since the Ni are bounded oper- 
ators, the condition TrFiV = N defines a closed linear 
subspace of £i{H). The intersection of this subspace 
with £^^{0, M] is therefore compact. £ is lower semicon- 
tinuous by fLemma 13. ip . Since a lower semicontinuous 
function on a compact space is guaranteed to have a min- 
imizer, the theorem is proven. □ 

A cg-density po is cg-EV-representable if and only if 
some r e S{po) is a ground state of some v ^ V. Since 
Theorem l3 . 2l shows that S_{po) is nonempty, the condition 
for cg-EV-representability is existence of G V such that 
^l(Po) + / vpo dx < Fi,{p) + J vpdx, for all p. This con- 
dition rearranges to the statement that the hyperplane 
p -Fl(Po)+ / v{p~Po) dx lies below the graph of Fl and 
touches it at po- That is, t; is a tangent functional to Fl 
at Po and, since it is continuous, a subaradient^. [ill, [l^ 
at Po- Recall that the set of all subgradients at p [de- 
noted 9Fl(p)], is called the subdifferential of Fl at p, 
and FL is said to be subdifferentiable at p if dFi^(p) ^ 0. 
It is a basic fact (ET: Prop. 1.5.2, vT: Thm. 5.35) that 
a convex function on a finite-dimensional space is sub- 
differentiable everywhere on the relative interior of its 
effective domain. This observation completes the proof 
of the central result of this section: 

Theorem 3.3 (Ensemble w-representability) . Every 
p in is cg-EV-representable. 

Proof. See preceding discussion. □ 

Combining this theorem with the cg-Hohenberg-Kohn 
Theorem gives Gateaux differentiability [ill, Il3| | of the 
restriction of Fl to Xn- 



Theorem 3.4 (Differentiability of FlIa-^). Fl, as a 

function on , is Gateaux differentiable on X^^ . 

Proof. Given p G X^'^ , Thm. 13.31 savs that there is at 
least one vinV which represents p. On the other hand, 
Theorem 13.11 asserts that there is at most one. So, there 
is precisely one. Uniqueness of the subgradient at p to- 
gether with local boundedness of Fl shows (ET: Prop. 
1.5.3) that FlIa-jv is Gateaux differentiable at p. □ 

Note in passing that we have no comparable theorem 
for E{v). If V has distinct ground states which lie in 
different strata of 5, then E is certainly not Gateaux 
differentiable there. 



4. UNCONFINED SYSTEMS 

Of course, the box used in the previous section is a 
rather artificial element, so it is desirable to prove analo- 
gous results with all space accessible to the particles. We 
refer to this as the unconfined case, and in this section 
we will show cg-EV-representability of all normalized, 
strictly positive, cg-densities for the unconfined case. 

Similarly to before, is partitioned into a regular, 
countably infinite, array of cells Bi, B2, . ■ ., each of vol- 
ume fl. This time, the regularity serves the important 
purpose of ensuring uniform upper and lower bounds on 
the cell sizes. As in Section [31 S denotes the set of 
normalized mixed states, X^ the set of normalized cg- 
densities, and X^~^ those which are nowhere zero. Now, 
however, X is £^ , the space of sequences {pi,p2, ■ ■ ■) such 
that IIpII — J2\Pi\ < This space is a Banach space 
when equipped with the norm || p|| . The Lieb functional is 
defined as before, and has the same upper bound p^"^'^^ 
on X^^ , outside of which it is -\-oo. X'^ is contained 
in the closed subspace Xn ^ {p ^ £^ : Pi — -^Is 
however, the interior of X^ relative to Ajv is empty, not 

'^N ■ 

Naturally, we then consider potentials in the dual space 
to which is the space £°° of sequences {vi,V2,---) 
which are bounded: sup \vi\ < 00. However, we shall 
eventually find ourselves looking for potentials in the 
larger set 

V = {v = {vi,V2, . ■ .) ■ inivi > -00} , (9) 

of potentials which are required only to be bounded be- 
low. Note that the cg-Hohenberg-Kohn Theorem 13.11 
works for potentials in V. 

From the experience of the box case, one might expect 
the work here to split into two principal parts: proof of 
the existence of Levy search solutions, by demonstration 
that Fl is lower semicontinuous and that the search is 
taking place in a compact space, followed by proof of 
subdifferentiability using general results from the duality 
theory of convex analysis. The first expectation is cor- 
rect. Lemma |3. II goes through unaltered, but Lemma 13.21 
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will require a more sophisticated counterpart. The run- 
up to Theorem 14.21 uses some ideas from Lieb(6i] (par- 
ticularly Thm. 4.4 there), but the proof given here is 
self-contained. However, what convex analysis gives us 
in this case is limited to showing that -Fl(p) and E{v) 
are a Legendre-Fenchel transform pair (up to signs). Es- 
tablishing existence of a potential v which has ground 
state with density p G takes an ad hoc argument 

following the pattern of the proof by Chayes, Chayes and 
RuskaiiSj] of EV-representability on a lattice. The basic 
reason for this is that some of the potentials needed do 
not lie in the dual space to In fact, they are not even 
linear functionals, because they take the value -l-oo on 
some densities. 

The immediate goal is the general-purpose Thm. 14.11 
below, which deals with (fine-grained) densities. Call a 
set A of densities tight if, given e > 0, there exists a 
bounded set such that no p £ A puts more than a 
total of e particles outside A^. Rephrasing, A is tight if, 
given a tolerance, there is a phantom box which contains 
all the densities to that tolerance. The set of mixed states 
whose densities belong to A is denoted by 5(A). The 
theorem says, roughly, that a set of phantom boxes is 
enough to guarantee that the low energy states among 
S{A) are an almost finite-dimensional set, which is one 
of the important services the real box of the previous 
section performed. 

Some notation that will be used is gathered here: For 
A a bounded region, denote the Hilbert space projection 
corresponding to all particles being in A by Pa- ^^o + 1 is 
bounded below by 1, and therefore has a unique positive 
square root denoted by h which is also bounded below 
by 1. is a bounded operator. 

The next lemma is key. It is followed by three auxil- 
iary results, the proofs of which can be skipped without 
detriment to understanding of the main results. 

Lemma 4.1. For A a bounded region of space, 
P\£^^(0, M]P\ is compact. 

Proof. Let (PAr„PA)5^Li be a bounded sequence in 
P\£~^{0, M]Pa. We will show that there is a norm con- 
vergent subsequence. Since r„ G £~^{0,M], 7„ = hTnh 
is also a bounded sequence in CfCH). Now, Ci{TC) is 
the topological dual to the space of compact operators 
on 7i, and any bounded sequence in Ci{Tt), in partic- 
ular (7„), has a weak-* convergent subsequence. [This 
is a slightly subtle point. It requires both the Banach- 
Alaoglu theorem and the fact that since H is separable, 
the space of compact operators on Ti. is also separable. 
See, e.g., Thm. 3.17 of Rudin[14].] By relabeling, write 
the weak-* convergent subsequence again as 7„, so for 
any compact operator A, Tr7„A — > Tr7^, where 7 is the 
limit. But, the Hamiltonian Hq is locally compact [15], 
i.e., PaH^^ and h^^P\ are compact operators (Lemma 
14. 4[) . Since the product of a bounded operator and a 
compact operator is compact, for any bounded oper- 
ator B, TtBPaT^Pa = Tr {(BPA/i-')7n(/i-'PA)} ^ 
TiBPAh-^jh-^PA. Defining F := h'^jh'^, this says 



PaF^Pa^-PaFPa with respect to the weak topology (not 
weak-*!). Now apply Lemma [4.31 to deduce PaF^Pa — * 
PaFPa. So, Pa£^^{0, M]Pa is relatively compact; but it 
is also closed, by Lemma [XT] □ 

Lemma 4.2. For F G Ci{TL) and P an orthogonal pro- 
jection on n, if ||PFP||t^ > (1 - e79)||F||t^, and e < 3, 
then ||F - PFP||t^ < e\\T\\tr. 

Proof. Without loss, assume F is normalized, so F ~ 
J2 Cm\4>m){'t>m\, whcrc Cm > 0, X] Cm = 1 and (/>m are unit 
vectors. Write F = PFP + PTP^ + P-^FP + P-^FP-^, 
with P^ = 1 — P, and bound the last three terms in trace 
norm as follows. 

IjP^FP^IIf,. = TrP^FP-^ = 1 - TrPFP < e^/9, 

where the final inequality is by hypothesis. Using the 
triangle inequality followed by the Cauchy-Schwarz in- 
equality and then the above bound on ||P-'^FP^ ||tr, we 
get ||PFP-L||,, = ||P^FP||t, < ||P0™|| < 

(Ec™||P^0™|P)'^'(Ec™|iP^™f)'^' < e/3. Gather- 
ing the pieces, and using e < 3, gives the result. □ 

Lemma 4.3. // (F„)^]^ is a sequence in CfCH), then 
F„ T if and only if F„ F, and in that case F G 

Proof. The implication from norm convergence to weak 
convergence is trivial, so suppose F„— >F weakly. The 
cone Ci{TL) C Ci{H) is norm closed and convex, hence 
weakly closed. Thus, any weak limit of (F„)5" is also in 
C^in). Since ||F„||f^ = Tr 1F„ ^ TrF = ||F||t^, the case 
F = is easy, and we may as well assume the F„ and F 
are normalized. 

By the triangle inequality, 

||F„ - T\\tr < ||F„ - PF„P||t, + I|PF„P - PFPIIt, 

+ ||PFP-F||t,. (10) 

Here, P is a finite-dimensional projection, which can be 
chosen (depending on F) to make the last term smaller 
than e'^/72. By weak convergence, ||PF„P - PTP\\tr < 
for n greater than some M. By Lemma 14. 2[ 
n > M then implies ||F„ - PF„P||fr < e/2, since 
||F„-PF„P||t, < ||PF„P-PFP||t, + ||PFP-F||t,. 
Now plugging into Eq. pH]) completes the proof. □ 

The next lemma is not newp^. but proofs in the lit- 
erature are sketchy. 

Lemma 4.4. Hq + 1 is locally compact. That is, if A 
is a bounded region of M.^^ , and Qa denotes projection 
onto A, then Qa{Ho + 1)^^ is a compact operator. 

Proof. First, note that (1 + T + Ke)"^ = (1 + T)''^ - 
{I + Py'^Veeil + T + VeeY^ ■ Siucc the interaction is 
relatively bounded with respect to the kinetic energy, i.e., 
V"ee(l -\-T + Vee)~^ is bouudcd, it suffices to prove this 
result for Hq — T. 
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1 + T = (1 - V^) = h"^, in appropriate units. The 
inverse of /i^™ has integral kernel 



(1 + A:2)m ^2TTy 



Gm{x, y) is bounded by c\x — ypm-sw \x — y\ — * 0, as 
follows by scaling, and for large |a; — ?/|, Gm{x, y) falls off 
exponentially. 

Thus, for m large enough, the Hilbert-Schmidt norm of 
QK{x)Gm{,x,y), given by j^^f^\Gm{,x,y)\'^ dx dy is finite, 
so Qa/i"^*" is Hilbert-Schmidt and a fortiori compact. 
Thus, h?"^ is locally compact. 

Now, if A > 1 is locally compact, then so is A^/"^. 
To see this, let u„ be a sequence of vectors converging 
weakly to zero. To show that A^^^^Qun 0, note 



that p-i/2gt 



{Qun\A ^Qun). By assumption. 



A-^Qun 0, so also \\A-^/^Qun\\'^ 0. Thus, local 
compactness of positive operators is preserved by taking 
square roots. Starting from the locally compact large 
power , repeated application shows h also to be lo- 
cally compact. □ 

Theorem 4.1. If A is a tight set of densities, and M < 
oo, then S{A) n £^^{0, M] is relatively compact. 

Proof. Choose sets Ai/„ corresponding to A as described 
before LemmaO Take T G 5(A) nf-^ (0, M], and note 
that r — Pai/^FPaj^^ is the part of the state which puts 
at least one particle outside Ax/„, so that its norm is 
less than 1/n: ||r — PAi/„rPAi/„ litr < l/n. This shows, 
via Lemma liTTl that S{A) f) £~^{0,M] is within l/n of 
the totally bounded set Pa^^^£-^{0, M]Pai/„- bmce n is 
arbitrary, S{A) n £~^(0, M] is itself relatively compact. 

□ 

With Thm. 14. II and Lemma [??T] in hand, it is now easy 
to prove the main theorems. 

Theorem 4.2 (Existence of Levy search solutions). 

For every p G , S_{p) ^ 0. That is, there is a mixed 
state r G S{p) such that Tt HqT = F^ip). 

Proof. The set of densities leading to cg-density p is a 
tight set. Apply Thm. O to see that S{p) = S{p) n 
£^^(0, Pl™"^] is compact. Since £ is lower semicontinu- 
ous fLemma l3.1l) . it attains a minimum on S{p). □ 

Theorem 4.3 (Lower semicontinuity of Fl)- Pl is 

weak lower semicontinuous on . 

Proof. It is required to show that Um ~ {p ■ Fl{p) < M} 
is weak closed for any M. Since Um is convex (be- 
cause Fl is), it suffices to show that Um is norm closed. 
So, suppose that the sequence (pn)J^Li is contained in 
Um- Choose r„ G 5(p„) and apply Thm. 14.11 with 
A — {pi,p2,...} to see that {r„} is contained in a 
compact set, hence contains a convergent subsequence: 
TriQ —>■ T- Since the map from states to cg-densities is 
norm continuous (each pi is a bounded operator), this 
means that F G S{p). By lower semicontinuity of £ 
(Lemma EH), £(r) < M. Thus, Fl(p) <M. □ 



The ground state energy for v € £°° is the infimum of 
the spectrum of F[{v), also given by 

E{v) = inf{FL(p) +vp:pe £^}. 

It can happen that there is no state which actually 
achieves the infimum. This definition actually makes 
sense for d G V. Also, v is bounded below in terms 
of E{v): 



Mv, > E{v)-Fi: 



(11) 



Since Fl is clearly convex, and weak lower semicon- 
tinuous by Thm. it follows (ET: Prop. 1.4.1, vT: 
cor. to lemma 6.12) that Fl (on l^) and E (on i"^) are 
a Legendre-Fenchel transform pair, up to signs. That is. 



Fl(p) = sup{£;(t;) -V p :ve 



(12) 



Again, as in the confined case, p is cg-EV-representable 
if and only if there is i; G V such that such that Fi^{pq) + 
Jvpodx < Fi,{p) + J vpdx, for all p. Note that this 
is not the same thing as saying that the supremum on 
the right-hand side of Eq. is attained. That would 
mean that the cg-density in question was a ground state 
cg-density of a potential in £°°; not all are so. 

It is this last observation that requires the method 
of proof of cg-EV-representability to differ substantially 
from that employed in Section [3] In a remarkable pa- 
per, Chayes, Chayes and Ruskai^] investigated the repre- 
sentability problem on a lattice and showed that all den- 
sities are EV-representable. Surprisingly, their method 
of proof can be adapted to the current situation with 
only minor modifications. The general idea is that the 
potential which represents p should be the one which 
maximizes E{v) — p v. Find a maximizing sequence d", 
and, using a diagonal construction and the countability 
of the set of cells, extract a subsequence which converges 
on each cell, to a candidate potential v. Cell- wise con- 
vergence is too vague to ensure much about H{v), so 
the second part of the proof establishes that H{v°') ac- 
tually converges to H{v) in strong resolvent sense, which 
is enough to show that it has a ground state in S{p). 

Specific results from Reed and Simon [17| or 
Weidmann[l8l| are indicated by RS or W. 

Theorem 4.4. If p £ '^n^ : then p is cg-EV- 
representable by a potential in V. 

Proof. A constant can always be added to v to make 
E{v) = 0, so 

FL(p) = sup{-t;-p : E{v) = 0}. 

Let v" be a maximizing sequence for the right-hand side. 
That is, 



plv^ldx 



)j°<0 



p\v''\dx ^ Fi^ip). 



(13) 



u°>0 
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Since E{v°') = 0, < somewhere. Define, for each a, 
a cg-density 



J 9{-v°')pdx 



Then, 



prna. > F^{p^) > [ 

J 



\dx > 



<o 



p\v 



dx. 



t)°<0 



So, J^Q^(j is bounded, uniformly in a, hence, 
by Eq. dill), /^^^^ is, too. Thus, / p|t)"| is 
bounded by some constant, and 



< 



c 



(14) 



Since each thus hes in a bounded interval, a diagonal 
construction will yield a subsequence which converges 
on each cell, though not in general with any uniformity. 
By relabeling, we can just write vf for this subsequence, 
so vf —>■ Vi for each i. This defines our candidate poten- 
tial V to have p as ground state cg-density. 

H{v) makes sense because Vi is everywhere bounded 
below by Eq. pT|) . and everywhere finite by Eq. (fT4)) . 
so V E V. Let V be the set of smooth pure states with 
density nonzero only in a finite number of cells; since 
the Coulomb repulsion is positive and locally L^, V is a 
common domain of essential self-adjointness for all H{v') 
with v' e V (RS: Thm. X.28, Hislop and Sigal: Thm. 
8.14). For (f) a unit vector in D, < (0|if(v")0) since 
E{v°') = 0, so 

< lim {(j)\H{v")(j)) = {(t>\H{v)(t>) < oo, 

a — >oo 

because the convergence of the potential is uniform on 
the support of </>. Thus, since 2? is a core, H{v) > 0. 
Further, \\[H{v) - H{v")](j)\\ = ||[t;" - v]<j)\\ 0, which 
shows convergence of H{v°') to H{v) on V. This suf- 
fices (RS:VIII.25 or W:Thm. 9.16) to establish strong 
resolvent convergence. 

Now, find Fq £ S_{p), so that 



f (Fo) = F^{p) 



lim 

a — ^oo 



-pw" dx. 



A Fq which satisfies the first equality exists by Theorem 
14. 2( and the second equality holds by Eq. (fTS]) . Since 
is a bounded function on x > 0, e^^'^""^ e"^*^'"-' in 
strong operator topology (RS:VIIL20, or W:Thm. 9.17). 
Then h\v) > imphes that 



1 > TrFoe"^'") = lim TrF 



oe 



-Hiv") 



> e"^(^°) lim exp 



— / pw" dx 



= 1, (15) 



by use of Jensen's inequality and the fact that {v") is a 
maximizing sequence for Fi^{p). Thus, TrFpe^^'-'"-' — 1, 
which, since H{v) > 0, shows that Fq is concentrated on 
the ground state manifold of H(v). □ 



5. INDEFINITE PARTICLE NUMBER 

Open systems [l^ are also of interest for density func- 
tional theory. In this case, the total particle number can 
fluctuate - the system is in equilibrium with a particle 
reservoir at some chemical potential, which we take to be 
zero without loss of generality. The theory of the previ- 
ous two sections extends in full to this case of indefinite 
particle number. We will explicitly discuss only the un- 
confined case. 

The relevant state space is the Fock space 



where T-Lm denotes the m-particle Hilbert space. The 
set of states which interest us is not the entirety of 
Ci'^{T), but the subset S of states which are diagonal 



in particle number, i.e., those of the form ^ 



m=0 "'m 

where F*^™-* is a normalized m-particle mixed state, and 
X] Om — 1- Below, we write F^™-' for the normalized m- 
particle part of any mixed state F. The following are nat- 
ural generalizations of previous definitions to the Fock- 
space setting. The subset of states in S which have den- 
sity p is denoted S{p), and the grand canonical density 
functional is defined 



F{p) :=inf{TrFiJo: T G 5(p)} , 



(16) 



where is simply the direct sum of the kinetic-plus- 
interaction Hamiltonians for all particle numbers. No- 
tice that the mixed states in Eq. (fT6|) are not restricted 
to two consequtive total particle numbers, in contrast 
to the pro cedure introduced by Perdew, Parr, Levy and 
Balduz|19|. The ground state energy for potential v is 
the infimum of the spectrum of E{v) and a state which 
realizes it is a ground state. As before, we can also write 

E{v) inf {F{p) + v ■ p : p e i^} , 

and a ground state is a state which realizes the infimum. 
Now, the additive constant is important. Decreasing it is 
likely to result in a ground state with more particles. For 
some potentials, even though bounded below, E{v) = 
—oo because the energy can be lowered indefinitely by 
adding particles. 

The proofs of Lemma 13.11 and of the Hohenberg-Kohn 
Theorem apply to the current situation with almost no 
essential changes. The density which is everywhere zero 
is the exception. It is certainly EV-representable, but 
not uniquely. 

We now sketch the required alterations to the proofs 
of Sec. m Consider generalizing Thm. l4.1l to a tight sets 
of densities A all with the same total particle number N . 
5(A) now refers to the corresponding set of A^-diagonal 
Fock-space mixed states. Choose e > 0. If TrFiVtot — N, 
then clearly J2m>N/t '^"i < £■ Similarly, if £(T) < M, 
then X)me/ "™ < where 

/ = {m e N : £:(F(™)) > M/e}. 
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According to Thm. 14. li the set Am of m-particle mixed 
states with energy less than or equal to M/e and density 
less than or equal to some density in A is compact. Since 
5(A) is within 2e of the convex hull of Um<jv/e^m, and 
e is arbitrary, S{A) r\£~^{0, M] is itself totally bounded, 
i.e., relatively compact. 

With Lemma [XT] and Thm. I4.1l thus extended to Fock 
space, the proofs that 5(p) ^ and that F is lower 
semicontinuous proceed just as in Thms. 14.21 and 14.31 
Thus, we can write, as before, 

Fi^{p)^sup{E{v)~vp:v (17) 

What is different about the situation now is that the 
effective domain of Fl is all cg-densities: positivity is 
still required, but normalization is not. 

The final stage is to adapt the proof of Thm. 14.41 to 
show that any cg-density, everywhere strictly positive, is 
Fock-space cg-EV-representable. The trick of adjusting 
the additive constant in the potential to obtain E{v) ~ 
had a significant part in the proof of Thm. 14.41 where it 
was both a labor-saving device and a means of assuring 
that the potential did not drift off to infinity. In the Fock 
space setting, it is neither possible nor needed, and we 
proceed as follows. 

Given p, find a sequence of potentials such that 

E{v^)-v"-p /Fip). 

Since F{p) > 0, for a large enough 

E{v^) >v" p> vZ„,N, (18) 

where — infivf. But also, by considering a state 

with any number M of particles all located in a cell where 
the potential is favorable, E{v°') < c' +VminM . Com- 
bining these inequalities, 

c'JVP 

Choosing M to be the next integer larger than N, for 
example, this shows that v!^^^ is bounded below inde- 
pendently of a for a large enough. And, inserting this 
bound back into inequality (|18p yields a lower bound for 
the E{v°') as well. Thus, by a diagonal argument as be- 
fore, a subsequence (written without loss) exists for 
which v" converges to Vi for each i, and E{v'^) to E^o as 
a ^ oo. 

The rest of the proof proceeds largely as before, with 
a few Eoo inserted. The appropriate domain of essential 
self-adjointness V is now the set of states with density 
supported in a finite number of cells, and also with com- 
ponents in only a finite number of Tim- The inequality 
showing positivity of H{v) becomes 

Eoo < {<t>\H{v)(f) < oo. 



Since {E{v^), E{v'^), . . . , Eoc} is bounded below, 
{H{v^), H{v^), H{v)} is also, so e'^^"") e-"'^'"^ 
in strong operator topology. With Tq e S{p), as before, 
the counterpart of Eq. is 

g-fi^ > TrFoe--^^^) = . . . = e'^^. 
6. CONCLUSION 

At the simplest level, the results presented in this 
paper merely increase our knowledge about the V- 
representability problem. Within every coarse-grained 
equivalence class of densities, there is at least one which 
is EV-representable. In addition, there is one such which 
is representable by a unique potential which is constant 
on cells, and bounded below, though not necessarily (in 
the unconfined case) bounded above. Extensions of these 
results to spin density functional theory and multicom- 
ponent systems is straightforward. 

Valuable as those observations may be, that way of 
putting it misses the bigger idea that a coarse-grained 
point of view may be a generally superior way to look 
at the foundations of density functional theory. From 
that perspective, what has been achieved is support for 
the idea in that EV-representability is shown to not be 
a problem at all in a coarse-grained framework. Little is 
lost by this change, since the coarse-graining scale can be 
chosen as the size of an atomic nucleus, for example, and 
most knowledge we may have of the fine-grained level can 
be incorporated relatively easily. 

Establishing EV-representability is the first step to- 
ward controlling the local structure of the Lieb func- 
tional. Because, if Fl is diffcrcntiable at p, in any 
sense, the derivative is given by the representing poten- 
tial. With regard to this question, the coarse-grained 
theory has taken one step. We have shown that Fl is 
differentiable in the confined case, by taking advantage 
of the finite-dimensionality. The most immediate out- 
standing problem at this point is to demonstrate some 
sort of differentiability in the unconfined situation. Gen- 
eral principles avail us nothing because X'^ has empty 
interior relative to X]\[. (The fine-grained theory also has 
this problem [2^. A widespread view that Fl is Gateaux 
differentiable at EV-representable densities is in error.) 
Hopefully, it will be possible to show that if 77^ 7^ for 
only finitely many i and ^rji — 0, then Fl(p + st]) is 
differentiable as a function of s at s = 0. This property is 
weaker than Gateaux differentiability, but probably not 
in any important way. 
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